TRANSFER MATRIX EIGENVALUES OF THE 
ANISOTROPIC MULTIPARAMETRIC U MODEL 



O 

o 
o 

(N 
(N 



o 



(N 

cn 
oo 
O 

o 
o 



Jon Links^ and Angela Foerster^ 



^ Centre for Mathematical Physics 

Department of Mathematics 
The University of Queensland, 4072 
Australia 
e-mail jrl@maths.uq. edu. au 

"^Instituto de Fisica da UFRGS 
Avenida Bento Gongalves 9500 

Porto Alegre, RS - Brazil 

e-mail angela@if.ufrgs.br 

Abstract 

A multiparametric extension of the anisotropic U model is discussed which maintains integra- 
biUty. The i?-matrix solving the Yang-Baxter equation is obtained through a twisting construction 
applied to the underlying Uq{sl{2\l)) superalgebraic structure which introduces the additional free 
parameters that arise in the model. Three forms of Bethe ansatz solution for the transfer matrix 
eigenvalues are given which we show to be equivalent. 
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1. Introduction 



The quantum algebras, including the Z2-graded analogues known as quantum superalgebras, play 
a central role in the construction of solutions of the Yang-Baxter equation which in turn may be used 
to construct integrable one-dimensional quantum models. Quantum algebras as defined originally by 
Jimbo and Drinfeld 0, || arise as one-parameter deformations of the familiar Lie algebras in such a way 
that the resulting algebraic structure is that of a quasi-triangular Hopf algebra. The supersymmetric 
generalizations are defined in 0, |, || and the quasi-triangularity of these Hopf (super) algebras is 
discussed in 0. The importance of this class of algebras is the existence of a universal element, known 
as the i?-matrix, which gives a solution for the Yang-Baxter equation. For each given solution of the 
Yang-Baxter equation, there is a well known procedure called the Quantum Inverse Scattering Method 
(QISM) by which an integrable one dimensional quantum system is obtained. One of the key steps 
in the QISM is the construction of the transfer matrix which yields a family of mutually commuting 
operators (including the Hamiltonian) which are taken as constants of the motion for the system. 
Diagonalization of the transfer matrix is the main objective in the solution of the model and from this 
result many properties such as the ground state structure and elementary excitations can be deduced. 

Through the construction of twisting by elements satisfying the two co-cycle condition it became 
apparent that new quasi-triangular Hopf algebras could be manufactured from existing ones , and if 
the twisting procedure introduced additional free parameters then these would appear in the Hopf alge- 
bra structure |Q. Viewed in another perspective, there exists classes of quasi-triangular Hopf algebras 
in which each member of the class is related to all others via twisting. In terms of representation theory, 
this means that representations of all quasi-triangular Hopf algebras in a given class are equivalent. 

On the other hand, in applications to the area of constructing integrable quantum chain models 
some subtleties emerge. In a recent work, it was shown that in the case of Reshetikhin twists [^], which 
will be defined precisely later, the periodic multiparametric chains can be mapped to the standard ones 
with the inclusion of a generalized twisted boundary condition [jl0| . In that work, two physical examples 
of the anisotropic t — J and U models of correlated electrons (both derived through representations of 
the quantum superalgebra Uq{sl{2\l))) were constructed and also the Bethe ansatz equations obtained. 
From these solutions, it is evident that the additional parameters introduced by twisting cannot be 
transformed away and thus do impact on the physics that these models describe. 

Subsequent to the work of Reshetikhin 0, Engeldinger and Kempf [|ll| gave a more general pre- 
scription for the twisting element (or twistor) by relaxing the triangularity property. This latter 
construction opens the possibility for even more parameters to be incorporated into the model in an 
integrable fashion. When applied to the fundamental representation of the Uq{gl{m\n)) series, the 
method of Engeldinger and Kempf reproduces the construction of Reshetikhin. For this reason, we 
cannot use these more general techniques to obtain a more general extension of the anisotropic t — J 
model than that already obtained in |]10[. However, when applied to "higher spin" representations, such 
as the case to be considered here, the differences between the constructions begin to emerge. We will 
show that these additional parameters describe local basis transformations for the local Hamiltonians. 
This result will be illustrated for the anisotropic U model which will be discussed in some detail. Being 
a fermionic model, it is necessary to derive the model in a supersymmetric formulation of the QISM. 

Unlike the case of integrable models based on non-graded algebras, there exist many examples of 
models with an underlying superalgebraic structure which admit more than one Bethe ansatz solution 
I2|, O, 111, Il5|, |T6|, O, IT8[ . It is generally accepted that the non- uniqueness of solution stems from the 



fact that the definition for a Lie superalgebra in terms of a system of simple roots is not unique 



For the case of the supersymmetric t — J model two forms of solution were known long ago from the 
works of Lai pO| and Sutherland [^. With the advent of the algebraic form of the Bethe ansatz a 
third form was discovered in [^, |TB[. Moreover, Essler and Korepin showed further in []T2| through an 
analytical argument that the three forms of solution for the supersymmetric t — J were equivalent. 
For the isotropic U model, as a result of sharing the same supersymmetry algebra (viz. s/(2|l)) 



as the supersymmetric t — J model, there are also three forms of the Bethe ansatz solution. For the 



transfer matrix eigenvalues, two forms were obtained by Pfannmiiller and Frahm [|T5| using standard 
methods. The third form was eventually discovered also by Pfannmiiller and Frahm [|16] but for this 
case one has to resort to more sophisticated Bethe ansatz techniques which were developed by Abad 
and Rios 



22 



Recall that the anisotropic U model was first introduced and solved via the co-ordinate Bethe 
ansatz in [^] as a generalization of both the Bariev model [24| and the (isotropic) U model given 
It was subsequently shown that the model is also obtained from an i?-matrix solution of 



m 
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the Yang-Baxter equation obtained through the one-parameter family of four dimensional (minimal 
typical) representations of the quantum superalgebra f/q(s/(2|l)) p6[ (see also [^, 



11,11,13). By 



applying the twisting construction to this solution of the Yang-Baxter equation we will give the explicit 
form for the resulting i?-matrix and in turn derive a multiparametric generalization of the anisotropic 
U model. 

One of the main objectives of this paper is to present all three forms of the Bethe ansatz solution for 
the multiparametric anisotropic U model including the Bethe ansatz equations and explicit eigenvalue 
expressions for the transfer matrix. Since the model is not based on the fundamental representation 
of its supersymmetry algebra, we will adopt a generalization of the procedure used in the case of 
the solution of the integrable spin 1 chain (Fateev-Zamolodchikov model) [^1], |3^, In the model 
considered here this involves determining the L-operator which acts on the mixed tensor space of the 
fundamental module with the four dimensional module which describe the local quantum states. A 

55] for classes of f/g(s/(2|l)) models 



similar procedure has recently been undertaken by Gruneberg |3J, 
which include the usual anisotropic U model. The necessary step for the multiparametric case is the 
determination of the appropriate multiparametric L-operator. In principle, it may be possible to work 
directly with the i?-matrix for the Bethe ansatz solution in an analogous way to that used by Ramos 



and Martins in the isotropic case ||3^ or even the analytic Bethe ansatz approach developed by Tsuboi 
for quantum superalgebras However, we will not consider these options in this work. 

Once the three forms of the Bethe ansatz solution have been obtained, we will proceed to argue 
that the transfer matrix eigenvalues are equivalent. Our approach is somewhat different than that used 
by Essler and Korepin for the t — J model [1^. We use the result that the Bethe ansatz solution for 



a Uq{sl{l\l)) (free fermion) model can be approached in two different ways and then proceed to show 
that the three solutions for the U model have their origin in the differing forms of these f/g(s/(l|l)) 
eigenvalue expressions. 

One of the sequent results of the Bethe analysis for the transfer matrix eigenvalues is that the 
eigenvalues of the quantum transfer matrix are obtained with little additional effort. The quantum 
transfer matrix approach, which has primarily been developed by Kliimper and collaborators I^S], ^ 



40| , has proved to be a powerful method to determine thermodynamic properties of an integrable model 
at finite temperature. In the final section of this work we will present the eigenvalues of the quantum 
transfer matrix for the anisotropic multiparametric U model. 



2. Quantum Inverse Scattering Method 



We begin by reviewing the fundamental features of the QISM. Let R{u) G End V he a solution 
of the Yang-Baxter equation 

Ruiu - v)RMR23{v) = R2z{v)RMRi2{u - v). (1) 

For full generality, we consider the cases also when V denotes a Z2-graded vector space. In such 
instances it is necessary to impose the following rule for the tensor product multiplication of matrices: 

{a®h){c®d) = {-lf^^''^ac®hd (2) 



for matrices a, b, c, d of homogeneous degree. The symbol [a] G Z2 denotes the degree of the matrix 
a. The monodromy matrix is defined 

T{u) = Rql{u)Ro{l-i){u)....Roi{u) 

from which the transfer matrix is given by 

t(M) = stroT(M). (3) 

A consequence of (|l]) is 

[t{u),t{y)]=Q, \/u,veC. 

Above stro denotes the supertrace taken over the auxihary space which is labeled by 0. 

In the usual manner the Hamiltonian associated with the transfer matrix is defined by the relation 



u=0 



H= t-\u)-^t{u] 
du 

Assuming regularity of the i?-matrix; i.e. 

i?(0) = P 

where P is the (Z2-graded) permutation operator, yields 

L-l 

H = ^ + hLi 

i=l 

where the local two site Hamiltonians are given by 

h = ^PR{u) 

Another observable operator which is readily obtained from the QISM is the momentum operator 
which is given by 

j9 = ilnt(0) (4) 
It is more useful for our purposes to work with the exponentiated form 

T = exp(— ip) 

= m 

= PlL----Pl?,Pl2 

which we call the translation operator. It satisfies the relations 

Th^i^i^T ^ = /;,(j+i)(j+2), ThiiT'^ = hi2. 

Clearly 

[T, H] = 

which reflects translational invariance of the periodic model. 

By using an available construction for obtaining multiparametric quantum algebras, it is straight- 
forward to obtain the associated multiparametric quantum spin chain. Below we will describe this con- 
struction and show that in each case we can effectively map the additional parameters to a "generalized 
boundary condition". 



3. Reshetikhin Twists 



Let {A, A, R) denote a quasitriangular Hopf (super) algebra where A and R denote the co-product 
and i?-matrix respectively. Suppose that there exists an element F & A ^ A such that 

(A®/)(F) = Fi3F23, 
(/®A)(F) = Fi3Fi2, 

F12-F13F23 = F23F13F12. (5) 

Then {A, A^, R^) is also a quasitriangular Hopf (super) algebra with co-product and i?-matrix re- 
spectively given by 

A^ = F,2AF{,\ R^ = F2iRF,-,\ (6) 

Throughout we refer to F as a twistor. 

The result stated above is a little more general than that originally proposed by Reshetikhin and 
is due to Engeldinger and Kempf |n[]. In the original work ||^ Reshetikhin imposed the additional 
constraint (triangularity property) 

F12F21 = / O / 

and in the case that {A, A, R) is an affine quantum (super) algebra Reshetikhin gave the example that 
F can be chosen to be 

F = exp ^ (Hi Hj - Hj ® Hi) (7) 

i<j 

where {Hi} is a basis for the Cartan subalgebra of the affine quantum (super) algebra and the (pij, i < j 
are arbitrary complex parameters. However, following the construction of Engeldinger and Kempf it 
is possible to choose 

F = expJ2(Hi®Hj)(Pij (8) 

which obviously gives a twistor dependent on more free parameters. Note that it is also possible to 
extend the Cartan subalgebra by an additional central extension (not the usual central charge) Hq 
which will act as a scalar multiple of the identity operator in any irreducible representation. 

Suppose that vr is a loop representation of the affine quantum superalgebra. We let R{u), R^{u) 
be the (super)matrix representatives of R and R^ respectively, which both satisfy the Yang-Baxter 
equation eq. (|lD. As R{0) = P then R^{0) = P result of d^). We may construct the transfer 
matrix 

t^{u) = stro(7r«(^+i)(J®Af)0 

= stro(<(w)i?5^_y(«). ...<(«)) (9) 

where Af is defined recursively through 

A£ = (J®/....® A^) A£„i 

= (A^®/....®/) A£_i. (10) 

Again the subscripts and 1,2,...,L denote the auxiliary and quantum spaces respectively and stro is 
the supertrace over the zeroth space. From the Yang-Baxter equation it follows that the multipara- 
metric transfer matrices t^{u) form a commuting family. The associated multiparametric spin chain 
Hamiltonian is given by 



u=0 



L-1 



+ (11) 



with 

= ^^^^ v^., 
= FhP-^. (12) 



The above construction allows us a means to incorporate arbitrary parameters (through the ipij) into 
the Hamiltonian without corrupting integrability. We will refer to such extra variables in the model 
as gauge parameters. 

Through use of (^ we may alternatively write 

t^{u) = stro (7r®(^+i) (/ ® Jl) [(/ ® A,.)(Fioi?oii^io)] (/ ® JlV) 

with 

Gi = FiLFi(^L~i)----Fi{i+i)- (13) 

We now define a new transfer matrix 

t{u) = JiH^{u)Jl 

= stro(7r^(^+i)(/®Ai)(Fioi?oii^oi')) (14) 

where we have employed the convention to let F denote both the algebraic object and its (super)matrix 
representative. Through further use of (|^) we may show that 

t{u) = stro {FwF2o....Fl^RMRo(l-i){u)....R^i{u)F^^\...F^I) . 

In this setting the operator T assumes the form 

T = Fi2^ F^^^ ...F^j^ F2iFsi...FliPil...PisPi2 

and the associated Hamiltonian is given by 

H = t-\u)^t{u) 
du 



u=0 



L-l 



hi,i+i + hb. (15) 



i=l 



where 

h = F^L.i)L ■■■F^lF^L_iy..FLihLiFj;l... F^j^_^^^FiL ....F(^L-l)L- 

It is thus apparent that the matrix Ji transforms the multiparametric Hamiltonian to one where the 
parameters appear only in a generalized boundary term h^. Although the boundary term is transformed 
by a global operator, there are some important properties which should be noted that show that the 
generalized boundary interaction behaves as a two site operator. The first is that 

[hb, = for z 7^ 1, L - 1. (16) 

Thus if we think of the local Hamiltonians as observables then we can still independently measure the 
boundary two site energies and those within the bulk. Also translational invariance is maintained; viz. 

[T, H] = 

and more importantly 

ThbT~^ = hi2, T/i(i_i)iT"^ = hb 



Th 



T- 



h 



(i+l)(i+2), 



for i ^ L — 1. 



(17) 



Consequently, such a model can still be interpreted as describing a closed chain system. This situation 
bears close similarity with the closed quantum superalgebra invariant models f^ . 



4. Jacobs-Cornwell Twists 



More recently, a new type of twisting 2-cocycle has been introduced by Jacobs and Cornwell in 
their work on relating non-standard quantum algebras to standard ones. In notation as above, suppose 
there exists F & A® A which satisfies the following relations 

F12F23 = F23F12 
(A®J)F = F23F13 
(/®A)F = F12F13. 

Then {A, A^, R^) is also a quasi-triangular Hopf algebra with A^, given by eq. (||) above. 

As in the case of the Reshetikhin twists above, one may use the Jacobs-Cornwell twists to construct 
multiparametric chains which can also be mapped to system where the extra parameters occur only in 
a generalized boundary interaction. For these models the transformation takes the form 

Jl = Gl-1----G2Gi 

Gi = Fj(j+i)Fj(j+2)....-Fji,. 
The translation operator under this mapping becomes 

F = F~j^F~^l_^y...F{2^FLi....F2iPiL----Pl3Pl2 

while the generalized boundary term hf, in the Hamiltonian is 

h = F[j^...F-\-^^^FLi...FL(L_i)hLiF£^\^-^y..F£^F(^L_i)L...FiL. 

The properties (0,|1^) also hold for the Jacobs-Cornwell twists. 

Finally, it is important to mention that the twistor eq. @ qualifies as both a Jacobs-Cornwell twist 
as well as a Reshetikhin twist. In general however, these two classes of twistors are inequivalent (see 

5. Symmetric twists 



As mentioned earlier, the antisymmetrization condition eq. (|^) originally imposed by Reshetikhin 
can be relaxed and we can consider the more general Cartan subalgebra twists of the form eq. 1^ 
which introduce more free parameters into the Hamiltonian. We will show here however that these 
additional parameters simply describe local basis transformations and as such to not have a bearing 
on the spectrum of the model. 

We begin with the observation that all the twistors of the type eq. (^ close to form a commutative 
group. Moreover, each twistor can be expressed as a product of a symmetric and antisymmetric twistor 
through 

where 

F' = (Fl2.F2l)l/^ F'^ = {F,2.F,,'f'\ 



The square root in the above expression is well defined given that the twistors of the form eq. (H) are 
defined in terms of the exponential of an operator. Given two Cartan elements k, I, let us consider a 
symmetric twist given simply by 

F = exp{k(^l + Wk). (18) 
It is a straightforward exercise to show that 

F = exp{A{k.l)).Ui.U2 

where U = exp{—k.l) and consequently for a twistor of the type eq. ([T8| ) we simply have 

= UiU2....UlHU{^U^\...U^\ (19) 

The fact that all symmetric twists arise as products of twistors of the form eq. ([I8| ) allows us to conclude 
that the relation eq. (|T^) is true for any symmetric twistor. As a result, we do not expect that the 
Engeldinger-Kempf form of twistor eq. (|^) will not introduce more dynamical gauge parameters than 
twistors of the Reshetikhin form eq. (0). 

6. The quantum transfer matrix method 



The algebraic Bethe ansatz method for the solution of the multiparametric anisotropic U model 
which will be employed below also has an important application in the quantum transfer matrix method 
PB| , |D| which we will briefiy describe here. 

As a result of the previous discussions, we can conclude that the QISM allows us to determine that 
the relation between the Hamiltonian and the transfer matrix is of the general form 

t{u) = T exp[uH + o{u% 

Let us also define 

t{u) = stro {Rio{u)R2o{u)....Rlo{u)) ■ 
A similar calculation shows that we may write 

t{u) = T-^ exp[uH + o{u'^)] 

and thus 

(ti-P/Cyti-P/C))''^' = exp[-pH + 0(1/ C)] 

where (5 = 1/ kT and throughout we will assume that C is even. Let us define two quantites U and U 
through the relations 

U = tr {t{-(3/C)t{-(3/C)f^\ 
U = sti{t{-(3/C)t{-l3/C)f'\ 

Both U and U can be thought of as partition functions for a classical two-dimensional lattice model on 
a torus which differ only in boundary conditions. In the thermodynamic limit £ ^ oo we can neglect 
such a difference in the boundary conditions and conclude that 

lim U = lim U = Z 



where Z is the partition function for the quantum system with Hamiltonian H (derived from the 
QISM) in the thermodynamic limit. 



At this point we define the quantum transfer matrix with inhomogeneity x to be 

Q{u) = stro {Rcoix - u)Roic-i){x + u)....Rti{x - u)Roi{x + u)) 

which forms a commuting family [Q{u), Q{v)] = by virtue of the fact that the Yang-Baxter equation 
(0) is expressible in the equivalent form 

Ri2{u - v)Rf{^{-u)Rfii-v) = Rfi{-v)Rtr{-u)R^2{u - v) 

where throughout stj refers to the supertransposition taken over the ith space. Note that we may write 

t{u) = stro • 

The quantum transfer matrix allows us to now express Z as 

Z = lim str (g(0)^) 

with X = —P/^- The thermodynamical properties are determined by the maximum eigenvalue of Q{u). 
Given that Q{u) forms a commuting family, the traditional Bethe ansatz methods can be applied for 
the diagonalization. 

7. Multiparametric anisotropic U model 



As an application of the above formalism, here we will introduce the multiparametric anisotropic 
U model derived from an i?-matrix obtained from the quantum superalgebra Uq{sl{2\l)). This super- 
algebra has simple generators {eo, /o, /io, ci, /i, hi} corresponding to the simple roots associated with 
the Cartan matrix 

' 1 
-1 2 



A 



It is worthwhile mentioning here that we work in the standard root system of one bosonic and one 
fermionic simple root. Another system of simple roots exists which we will not consider (see e.g. [|30| ). 
For a full definition of the algebra we refer to ||28|| . 

This algebra admits a non-trivial one parameter family of four-dimensional representations, which 
we label vr, given by 



vr(eo) 


= v[a 


el + ^J[a + l]e\ 


vr(/o) 


= V[a 


e{ + + \\e\ 


7r(/io) 


= a{el 


+ e^) + (a + l)(e^ + el) 


7r(ei) 


= -el 




<h) 


= -el 






= el- 


e^- 



(20) 

Above the indices of the elementary matrices e* carry the Z2-grading (1) = (4) = 0, (2) = (3) = 1 and 
we employ the notation 

— 



q-q 1 



Throughout we work under the assumption that q is generic. Associated with this representation there 
is a solution of the Yang-Baxter equation which is obtained by solving Jimbo's equations. The problem 



of obtaining this solution has been considered in |2^, ^ J30[| . Applying the construction we 
described earlier yields the following multiparametric solution of (|1|) 
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r>32 
itl4 








r>32 
^23 








n32 
^32 








r>32 









































n33 
^33 
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-0,24 

















n42 
^42 









































Rtl 








n43 
-^43 


















































r44 

-n.44 



where the non-zero entries are given by 

[n— o] 



Rll 

Dl3 

* pl2 

* d14 



-'^32 

* p22 
-"-22 



p24 

* p23 

* p23 
-K14 

p31 

* p33 
-"-33 

p31 
^13 

p34 
^43 

* p32 
p42 



42 

p44 
-fl,44 

* p42 
^24 

* p41 
^23 



[u+a] 

nU ["] 
y [u+a] 

2u 

" [tt+a] [n— a— 1] 
_„ „-l„n+l/2 [a]l/^[a+l]l/^M 
PIP4 y [u+a\[u-a-l\ 



-1 [«] 

2g-g2" + l-g-2a-l-g2t.+ l^g2^.-l 

(^qu+oL_q — u — oL^(^qu — a — l__q — u+<^-\-l'^ 

-r> r, ^~«+l/2l^m«+lF!M 
P2P4q [u+a][u-a-l] 

-1 M 

P2 [u+a] 
1 

_„ Ir.] 



[u+a] 

„u [g+i] 

y [u-a-1] 

^^-1/2 [a]^/^[a+l]^/^[u] 
P2 P3q \u+a][u-a-l] 

-1 -1 [u\ 
P2 P3P4 1^1^ 
[u+a+1] 
[u—a—1] 

_ „ jo+iL 

y [u-a-l] 

„-1„-1„-«-1/2M!^VH&1 
Pi PS y [u+a][u-a-l] 



Dl2 
-"-12 

pl4 

* pl3 
^31 

* pl4 
-^23 

p21 
^21 

* p23 
-'^23 

p21 
^12 

p24 
^42 

* p23 
-K4I 

* p32 
^32 



p34 
rt34 

* p32 
-'^23 

* p32 
p41 

D43 
^43 

p41 

* p43 

-n^34 

* p41 
^32 



Pi 



[u+a] 



[u][u-l\ 



P^P"^ [u+a][u-a-l\ 
„u [q] 
y [u+a] 

„ „-l„«-l/2 [a]^^^a+lY/M 

P2P3 y 

-1 [u] 
Pi [u+a] 

-1 -1 
Pi P2P3 P4 [u+a][u-a-l] 

y 

y [«-«-!] 

„»+l/2 [a]l/^[a+l]V^M 
Pi P4H [u+a][u-a-l] 
-1 ~1 [n]2 



2^-l_^2a+l_^-2Q,-l4_„-2ii+l_„-2u-l 



„ -«-i/2[an^>+i&] 

f'lf'3y [«+Q][«-a-l] 

-1 -1 M["-i] 

^'2 [u+a]\u-a-l] 
- 1 - 1 mI 

Pi P3 Piy;^] 

2u 

" [m+o] [u—a—1] 
y [m-q-1] 

„-l„-l„-»+l/2 [a]'/^[a+lV'^u] 
P2 P4 y [u+a][u-a-l] 



(21) 



and we adopt the notation 



y-y ^ 



We remind the reader that the above i?-matrix solves the Yang-Baxter equation eq. (|T]) subject to the 
rule eq. (0) which is a consequence of the f/q(s/(2|l)) superalgebraic structure underlying this solution. 



The significance of tlie * notation prefixing some matrix elements will be explained in a subsequent 
section. 

Using the above solution, we can employ the QISM to obtain an integrable Hamiltonian as discussed 
earher. Identifying the four basis states in the f/q(s/(2|l)) representation space V with the local 
electronic states through 

.;^ = |0), v' = \i), v' = n), v' = \n) (22) 
allows us to express the local Hamiltonians in the following form (with convenient normalization) 

-p-^p^%cy^,^,^^ {qpl[a + {q-'pf[a + 1] H-^) + h.c. 

+ H"^(PrV2^4T4iC(i+i)iC(i+i)| +h.c.) + [a]-\nii;ni^ + rii^Uii) 

+ - 1) + g-"-^(n(,+i)^ + - 1). (23) 

Above we have used standard notation; the operators c^, c denote fermi creation and annihilation 
operators and n measures occupation number. For the above operator to be hermitian, it is assumed 
that the parameters pj, j = 1,2,3,4 lie on the unit circle. This results from the fact that in the h.c. 
terms one has pj — > pj^. 

An immediate feature of this Hamiltonian is that it is not spin reflection invariant. However, it is 
invariant with respect to spin reflection coupled with the interchange of parameters 

q ^ pi ^ p2, P3 ^ Pi- (24) 

This invariance will manifest itself in the Bethe ansatz solutions determined later. 

Although the above local Hamiltonian depends on the four parameters pj, we can apply the unitary 
transformation 

^ 4t(P3P4)"^"'S 4^ cJ^(P3P4)^"'^ 

which yields the following local Hamiltonians 

= -P2" V3P4"'cItC(.+i)T {q'VpAc, + 1] [«]"') {QP3'P4[a + 1] [a]-') + h.c. 

-p^Wp4lCi^+l)l + ^"'^ {q-'p3P^> + + h.c. 

+ [a]-'^{p^^p^^cl^cl^C(i+i)iC(i+i)'i + h.c.) + [a]-'^{nii;nii + rii^rin) 

+g"+i(n,| + - 1) + g-"~'(n(,+i)^ + n(i+i); - 1). (25) 

It is clear that there are only three independent gauge parameters pi, p2, ^3^4^- This result is confirmed 
by the transfer matrix eigenvalue expression and Bethe ansatz equations, which will be derived below, 
as they exhibit dependency on only these three independent gauge parameters. In the context of our 
earlier discussion, the parameter pap^ enters into the Hamiltonian by means of a symmetric twistor 
and thus may be transformed away. 

A final comment here is that the above Hamiltonian is not the most general that can be obtained 
through this procedure. By taking the tensor product representation of Uq{sl{2\l)) with different 
values of the free parameter a in (^) a more general solution of the Yang-Baxter equation is obtained. 
For the gauge free case the explicit form is given in [p4| , Using the method proposed in a 
Hamiltonian can still be derived through this solution. 



8. Bet he ansatz solution 



Having derived the model, we now use the method of the algebraic Bethe ansatz in order to find 
the eigenvalues of the transfer matrix which in turn permits us to calculate a formula for the energy 
levels of the Hamiltonian. Since we are working in a "higher spin" representation of f/g(s/(2|l)), we 
follow the approach first proposed for the study of the integrable spin 1 chain |^T], 32, Rather than 
seek "creation operators" over the pseudo- vacuum from the action of R{u) (see e.g [^) we instead 
appeal to an algebraic structure derived from a lower rank object. In order to achieve this goal, we 
first introduce the matrices TZ{u) G End(l^ ® W) and L{u) G End(iy V^) where W denotes a three 
dimensional space and V is the four dimensional representation space of f/q(s/(2|l)) as before. We 
require these operators to satisfy the following forms of the Yang-Baxter equation 
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Explicitly, TZ{u) and L{u) satisfying the above relations take the form 
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The existence of these solutions is a consequence of the representation theory of the untwisted 
affine extension of Uq{sl{2\l)) (e.g. see |2^) where the three dimensional space corresponds to the 



fundamental module. Again, multiplication of tensor products in (P6| , P7| , P8D is to be undertaken with 
consideration to @. In order the simplify the task of the Bethe ansatz approach we work here with the 
solutions of (^,^,^) without the Z2-graded tensor rule. This is achieved by a change of sign by those 
matrix elements of TZiu), L{u), R{u) prefixed with a *. Again, we refer to for a justification of this 



approach. We comment that in the space W the Z2-grading of the basis states is (1) = (2) = 1, (3) = 0. 

Let us introduce two associative algebras Y and Z with elements {yj{u)}^j=i and {Zj{u)}fj^^ 
respectively and subject to the relations 

TZuiu - v)Yisiu)Y23 = Y23iv)Yr3{u)ni2iu-v), (31) 

Ruiu - v)Zi3{u)Z23{v) = Z23{v)Zi3{u)Ri2{u - v), (32) 

Luiu - v)Yi3{u)Z23{v) = Z23{v)Yi3{u)Li2{u - v). (33) 

where 

Y{u) = 5^e} ® Yj^, Z{u) = Y^e)® Zi{u). (34) 

The associativity of these algebras is guaranteed by the relations eq. (|l]) and eq. (^6]). By comparison 
with eq. (^ we see that the monodromy matrix provides a representation of the Z algebra acting on 
the module V®^ by the mapping 

vr {Z]{v))\ = (-l)«(0+0-)(0+W(fc)Tifc(^). (35) 
Moreover, the transfer matrix is expressible in terms of this representation by 

4 

1=1 

The phase factors present above are required since the Z is defined in terms of the non-graded L- 
operator. 

We also define an auxiliary monodromy matrix through 

U{u) = Lol{u)....Lq2{u)Lqi{u). 

Likewise, the auxiliary monodromy matrix U (u) gives a representation of the Y algebra through the 
action 

vr {Yjiu))'; = (-l)»{0+0-)(0+WWf/|(„) (37) 
while the existence of the solution for (l28| ) permits us to appeal to ( ^3]) for algebraic relations between 



elements of Y and Z. In the following we omit the symbol vr for ease of notation. 



We choose as a reference state (pseudo-vacuum) for the Bethe ansatz procedure the L-fold copy of 
the highest weight state in V; viz. 

$0 = (v'f^ (38) 
which itself is an eigenstate of the transfer matrix with the eigenvalue 
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In order to obtain more eigenstates of the transfer matrix we adopt the ansatz 



m 



(39) 



(40) 



where the F^^y are undetermined co-efficients and we have set 

S^^Hiu}) = Yt{u{)Yt{u,)....Yt{uM) A = 1,2. (41) 

The motivation for this choice of ansatz is standard in the sense that a routine calculation shows 

y/(«)$o = VjVz^3 (42) 

while is an eigenstate of Y^{u)^ i = 1, 2, 3. 

We now appeal to the algebraic relations given by equation ( |33[ ) in order to determine the constraints 
on the variables Ui needed to force equation ( ^OD to be an eigenstate. Of the numerous relations resulting 
from (P^D we need only consider the following. 
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In the above commutation relations two types of terms arise when the actions of the operators 
Zl[v), Zf[v) are evaluated on (0). These terms are those giving a vector proportional to which 
we call wanted terms while all other types are unwanted terms (u.t.). The operators in the above 
expressions contributing to the unwanted terms are all those of the form Zj{v)Y/'{u) with the necessity 
that i < j. Note that for any Yi^[u) term which arises where / 7^ 3 and A; 7^ 1, 2 we can appeal to 
the commutation relations (|3lD to commute it through to the reference state where it will act as 



a scalar. Repeated use of this procedure shows that all unwanted terms are ultimately expressible in 
the form 

s^''H{u'})z;{v)s^^'\{u"})^^ 

with i < j and v arbitrary and so can never give a vector proportional to . 

In terms of the fermionic representation (^2]), it is deduced that the reference state is the 
completely filled state with two electrons of opposite spin at each lattice site. Furthermore, it is not 
difficult to see that the operator Ll{v) annihilates a spin up electron from the system while L'^{v) 
annihilates a spin down electron. Using the condition that the states (0) have well defined electron 
number and spin (this can be interpreted as the conservation condition in the sense of Schultz p^), 
we have 
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where Uj^ and n_ are respectively the number of spin up and down electrons in the state and = 
n+ + n_. 

The situation for determining the action of ^K^) is a little more complicated. For these 

cases more than one wanted term may arise in the calculation, which forces us to employ a nested 
Bethe ansatz approach. In terms of the six-vertex solution of the Yang-Baxter equation with gauge 
parameter / 

([v + l] \ 

l[v] 
q'"" l-\v\ 
\ [^^ + 1] / 



r(t>, /) 



(43) 



we define the monodromy matrix 

t{v, {u}, I) = roi{v - ui, l)rQ2{v - M2, /)....roiv(w - u^v, /). 
It turns out that in this notation we may write 
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Completing the diagonalization of the transfer matrix is accomplished by diagonalizing the operator 

fiitl{v,{u},l) + fi2tl{v,{u},l) (44) 

which is the transfer matrix for the six- vertex model with gauge parameter /, inhomogeneities Ui 
and twisted boundary conditions. This is a standard calculation so we simply present the eigenvalue 
expression which reads 



2L~N 



M 



X{V) = /ii/*^ J] _ + 1] J] 
i=l j=l 



\v — w 



1] 



2L-N 



M 



V - Wj] 



i=i j=i 



W, + 1] 



Making the substitutions 



V - 
I - 

fJ'2 



> V — a/2, 
-^Pl'p^Ps'^P^ 



allows us to determine the contribution to the eigenvalue expression for the transfer matrix coming 
from the terms and Z^{v). At this point we would like to remark that the reference state used 

in the Bethe ansatz diagonalization of the six-vertex transfer matrix (^i]) corresponds to the state 

Lliu,)Lliu,)....LlM^' 

which is an 2L — electron state containing L spin down electrons and L — N spin up electrons. 
We can immediately deduce that M = L — n^. Collecting these results and making the substitutions 



Ui — i> Mj + 1, Wj 



Wj + 1/2, we may now give the eigenvalue expression 
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(45) 



for the transfer matrix (j^). Constraints on the parameters Ui, wj are given by the Bethe Ansatz 
equations. In principle, these may be derived from the requirement that all the unwanted terms cancel 
in the above calculation. A more practical way to obtain them is to use the fact that the eigenvalue 
expression ( ^S] ) is an analytic function of the variable v. Imposing this condition leads us to the set of 
equations 
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Alternatively, we can use the following reference state 

Ll{u^)Ll{u2)....LlM<l>' 
for the diagonalization of the operator (0). In this instance we obtain 
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for the eigenvalues of the transfer matrix. It is clear that (0) is obtained from (^5]) by the change 
of parameters (|2^ ) and spin reflection, a consequence of the invariance of the Hamiltonian p^ ) with 
respect to this operation. For this case the Bethe ansatz equations assume the form 

\[u, + a/2]J Al iu,-wj-l/2] 

2L-N r ^ L-n+ , 

ftftft ft 11 _ + 1/2) - n J- 

It is not clear whether the two eigenvalue expressions presented above, either individually or together, 
give the entire spectrum for the model. No highest weight theorem for the Bethe states is applicable 
here due to the fact that the imposition of periodic boundary conditions does not preserve invariance 
with respect to the action of the elements of Uq{sl{2\l)). In the q —>■ 1 limit the invariance is restored 
and one can show that the Bethe states are indeed highest weight states similar to the methods used 



for other s/(2|l) invariant models 0, |T^. However the combinatorial arguments given in those works 
do not extend in any obvious manner for this case and to our knowledge the completeness of the Bethe 
states for the U model remains an open question. 

Using the transfer matrix eigenvalue expressions ( P5| , |47| ) the energies of the Bethe states in both 
cases are obtained through 
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One can alternatively use a tensor product of lowest weight states as the pseudo-vacuum for the 
Bethe ansatz calculation; viz. 

This leads to a different form for the transfer matrix eigenvalues and Bethe ansatz equations. The 
procedure used is analogous to that above, and so we will omit the details of the calculation and just 
give the final results. The transfer matrix eigenvalues are 
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where the Bethe ansatz equations read 



Just as in the previous case there is a second expression for the eigenvalues, due to a different choice 
of reference state in the nesting, which is 
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For this expression, the Bethe ansatz equations are 
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One can see that the formulae A2(f,Mi,Wj) and K2{v,Ui,Wj) are related through the mapping (p^) 
and spin reflection as in the previous example. Moreover, there is a mapping from Ai(f,Mj,tfj) to 
A2(t', Mj, Wj) and Ai(t>, Wj, Wj) to A2(f , Uj, Wj) (and the corresponding Bethe ansatz equations) which is 
given by 

Pi P2^P3P^^, P2 PrVs V, P3 Pi, P4 ^ P3, a ^ -a - 1, n+ ^ L - n_, n_ ^ L - n+. (50) 



The existence of this relationship can be understood in terms of the dual representation of (pO]). Recall 
that for any representation of a quantum superalgebra vr the dual representation tt* is defined for each 
element a in the superalgebra through 



St 



Tx*{a) =7r(5(a)) 

where st denotes the supertransposition as before and S is the antipode. For the representation (pO|) 
under consideration here, the dual representation is equivalent (i.e. up to a basis transformation) to the 
change of variable a —>■ —a — 1. The change in the gauge parameters is also a result of the dualization 
of the Cartan elements in the definition of (|]). 

The transfer matrix eigenvalue expressions ( P5| , ^DD for this second form of the Bethe ansatz give 
the energies of the Bethe states through 



a+l _ -a-1 

E = — ^ A ^ — 

2 In g du 



N 



u=0 

[a + If 



^ [«, + «/2 + l/2][M,-a/2-l/2]' 



A third form of the Bethe ansatz also exists for this model, which was discovered by Pfanmiiller 
and Frahm [16| in the isotropic case (g = 1) with all gauge parameters equal to 1. The extension to 
the present model is achieved by a similar construction and so we will again omit the details and refer 
the interested reader to |T^ for the methodology. Again, there are two types of eigenvalue expression 
which are 



P2[v] 



n 



Pi[ui-v - a/2 + l/2\^ [wj-v + a/2\ 



n 



+ 



[v + a]J -Li [ui-v + a/2 + 1/2] J-A p2 K' - f - a/2] 

t — i .7 — i 



PiPsiv - a - 1] 
PiPslv + a\[v - a-1] 

L—n 



nf pm[ui 
\ r\ All ■ - 



v-a/2- 1/2] 
\ V Pi[ui-v + a/2 + 1/2] 



n 



n 

1=1 



n PiPsjuj - v-a/2 - 1/2] 



p^[Ui-V- a/2 + 1/2] 
[ui-v + a/2 + 1/2] 

[wj — V + a/2] 



Pslwj -v + a/2 + 1] 
P2PA[wj - v- a/2] 

p?J[Wj -v + a/2 + 1] 



n 



P2Pi[wj - v- a/2] 



^ Pi[ui- V + a/2 + 1/2] j^^p2[wj - V - a/2] 
subject to the Bethe ansatz equations 



— T. n+~L L—n+ 

Pi P3 Pa 



[uj + a/2 + 1/2] 
[ui - a/2 - 1/2] 



L n——L L—n 

P2P3 Pa 



for which the energy is given by 

L—n— 



K- -«/2] 
K- + a/2] 

[a + 1]2 



-Q [wj - Ui - 1/2] 
i= 

L — 71_ 



i |u,j - Ui + 1/2| 



n 



^ [ui + a/2 + l/2][ui- a/2- 1/2] [wj + a/2][wj - a/2] 



n+ 

E 



[uj - Wj - 1/2] 
\ [u,-w, + 1/2] 

[a][a + l] 



and 



L L-n+ 



Pl[V\ 

[v + a] 



n P2[ui - V - a/2 + 1/2] -p-r [wj - V + a/2] 
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+ 



\ [ui-v + a/2 + 1/2] Al [ujj - V- a/2] 

L L-n 



nf P2Pa[ 



Ui-v - a/2-1/2] 



P2Pa[v -a-1] J fj^ \ p3[ui - V + a/2 + 1/2] 

L L-nj, 



n 



PiPz[v] 



P2Pa[v + a][v - a - 1] 



n 



P2[u,- V-a/2 + 1/2] 
[u,-v + a/2 + 1/2] 



P4[wj - V + a/2 + 1] 
Pmlwj - v-a/2] 

P4[wj - V + a/2 + 1] 



(51) 



(52) 



n 



Pm[wj - v-a/2] 



P2PA[Ui -V- a/2 - 1/2] -|-r 
nJi/. — 7, -U /9 -U 1 /9l 11 



V + a/2] 



_^ P3[ui - V + a/2 + 1/2] jJ^Piiwj - V - a/2] 
in which case the Bethe ansatz equations are 



(53) 



-L L-n. n.-L ( [Ui + a / 2 + 1 / 2] 

P2 P3 Pa 



[ui - a/2 - 1/2] 
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L n+—L L— n.j 

P2P3 Pa 



[wj-a/2] 
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J 

L—nj^ 

n 

i=l 



[wj - Uj - 1/2] 
\ K- -Ui + 1/2] 

Ui - Wj - 1/2] 



\Ui 



w, + 1/2] 



(54) 



and the energy expression is 



^ [u, + a/2 + 1/2] [u, - a/2 - 1/2] ^ [w^ + a/2] K- - a/2] ' 

Again, the two eigenvalue expressions and associated Bethe ansatz equations are related through spin 
reflection coupled with the change of parameters (|2^) . A more striking feature is that ( PH , ^ , ^ , ^ ) 
are all invariant under (0) with an accompanying interchange of the sets of parameters {ui} and {vi}. 

9. Equivalence of the transfer matrix eigenvalues 



The above calculations show that there are three different forms for the Bethe ansatz solution of 
this model. Here we will argue that these forms are in fact equivalent. The method that we will use 
is inspired by the paper which hints at the argument we will give without providing many details. 
Our goal here is to make the argument more transparent. 

We begin with the following i?-matrix ^ which is associated with the quantum superalgebra 



Uq{sl{l\l)). The solution we present below is more general than that of 28 1 in that we have included 



a gauge parameter / which is necessary for our argument to succeed. The i?-matrix is 

j\v-{a + (3)/2-l] \ 

/[t; + (a-/3)/2] gV[« + l][/5 + l] 

^ -g-V[« + l][/^+l] l-\v + {P-a)/2] 

V [v + {a + p)/2 + l]) 

(55) 

which satisfies the coloured Yang-Baxter equation 

r1^{u - V, l)r^^{u, l)r(^3{v, I) = r^^{v, /)rJ7(M, l)r"^{u -v,l). 

We emphasize again that the tensor products in the above equation are to be evaluated according to 
the rule (H). In (|55[) we have adopted the convention that index 1 is odd and 2 is even. 



Using this solution we can construct the transfer matrix 

t{v, a, P,Ui,l, Xi,X2) = stro (^Xor^^{v - UiJ)...r^^{v - Up,l)r'^^p_^^^^{v,l)...r°^^_^p^{v,l)^ 
where X = diag(a;i, X2), x G C These matrices form a commuting family with the property 

[t{v, a, P, Ui, I, xi, X2), t{u, a, 7, Ui, I, Xi, X2)] = 0. 

For the diagonalization of the above transfer matrix by the Bethe ansatz there are two available 
approaches. We may begin with a reference state given by 

in which case the eigenvalues read 

/ P 

A {v,a,P,Ui,Wj,l,xi,X2) = lx2[v + {a + P)/2 + l]^Y\_ [v-Ui + a/2 + 1] 



i=l 

P \ M 



such that the parameters Wj satisfy the Bethe ansatz equations 

K- + 13/2 + 1/2] ,L /[Ui - + 1/2] 



X2 



x,l^W T~ ■ (57) 



w^- (3/2-1/2]) 11 [^^_^._i/2]- 



C^lN®(L+P) 



We can also use 
for the pseudo- vacuum which yields the eigenvalue expression 

A^(f , a, (3, Ui, Wj, I, Xi, X2) = — I Xi[v — {a + P) /2 — 1]^ Y\ [v — Ui — a/2 — 1] 



N 



i=l 

P \ L+P-M 



-.,r'|„ + (,-a)/2]^nr.|.-..-a./2|) U (68) 

such that the parameters Wj satisfy the Bethe ansatz equations 

W,+,3/2+l/2] y_ ,-A ik-W, + l/2] 
' |W,-,8/2-l/2|j - 11 K-W,-l/2| 



We conjecture that for each set of M parameters {wj} satisfying (|5^), there is a set of (L+P-M) 
parameters {wi} satisfying ( ^9]) which render (|56D and (|58|) equal. In the rational limit g — > 1 this 
argument can be made rigorous since in this instance there is an underlying s/(l|l) invariance for the 
transfer matrix. As mentioned before, one can prove a highest and lowest weight theorem respectively 
for the two sets of Bethe states and then use a combinatorial argument to claim that the multiplets 
generated by the Bethe states give the full space of states in each case, so there must be a one-to-one 
correspondence. Although this proof will fail for generic values of q, we still believe that the result 
holds true. 

Applying the conjecture, it then follows that 



As{v,Ui,Wj) 



L—ri- 



1 



[v + a]^ y [v -Ui- a/2 - 1/2] 



X < — A f , a — 1, a — 1, Mi, Wj,p2 



+ - —j-X \ V - 12, a, a - l,Ui,Wj, , 

[v -a-l\^ y P3 \ P-i J \ Pi 

L—n- ^ 



^— n 

[v + a]^ Al [v-Ui- a/2 - 1/2] 

i \+ ( 1 1 - f PiV^^'''^ /PiPsV^""^ f PiPsX 

X <-X^ \v,a-l,a-l,Ui,Wj,p2,\—j ' 1 J ' 1"^^/ 



[v]^ X+ 1 /O 1 - ~n. f P2P4Y'' f Pm 

\v - l/2,a,a - l,Ui,Wj,——,p^ V~P~/ 'V 



[v - a - 1]^ y P3 \ P3 J \ Pi 

Ai{v, Wj,Ui) (60) 

(61) 



Similarly, we find 

As{v,Ui,Wj) (62) 

L— n+ 



4- r^]L n 



+ a]L 11 [v-u,-a/2-l/2] 
X <( — A~ I f, a — 1, a — 1, Uj, Wj,pi, 



y —rX \v -12, a,a-l,Ui,Wj, ,P2 A 

[v-a-1]^ \ ' Pi \ Pi J \ P3 J 



L-n+ ^ 



W+W ^ [v-u,-a/2-l/2] 

\ (L— n+) / \ lij— n+i / 



\PlJ \ Ps J \ P3 



V P4 V P4 / \ P3 J 



[v — a — 1]^ 

Ai{v, Wj,Ui) (63) 

(64) 



although the result is immediate through use of (plf). In a similar fashion it can be shown that 

A3{v,Ui,Wj) = A2{v,Ui,Wj), A3{v,Ui,Wj) = A2{v,Ui,Wj) 

though the judicious mathematician would appeal to (|50|) to obtain the result. 



10. Quantum transfer matrix eigenvalues 



The Bethe ansatz approach for the diagonalization of the quantum transfer matrix proceeds in 
much the same way as the previous discussion. We return to the algebraic relations (|3T|j3^ , |33|) . A 
representation of the algebras Y and Z are obtained through 

niYiv)) = {L^l-ix-v)Loic-i)ix + v)....L%^ix-v)Loiix + v)) 

7r{Z{v)) = {R^l^(x-v)Roic-i){x + v)....R%^{x-v)Roi{x + v)) 

where the label refers to the auxiliary space (three dimensional for L{v) and four dimensional for 
R{v)) and the natural numbers label quantum spaces. The quantum transfer matrix is given by 

4 

Q{v) = Y.{-iY'^^^'^^'\{zi{v))'; 

i=l 

in complete analogy with (|36|). The principle difference is in the choice of the reference state. In order 
to adopt the ansatz (^) we need to find a reference state which satisfies the conditions (^). This is 
achieved with the choice 

which is an eigenstate of the quantum transfer matrix with eigenvalue 

[x + v - a][x-v][x-v \ _ f P3[x + v][x-v] \ 

,PiP2[x + V + a][x- V + a][x- V - a - 1]) \PiP2P4[x + v + a][x - v - a - 1]) 

Pi[x + v][x-v] f [x + v][x + v-l][x-v + a + l] ^ 



PmPslx + v + a][x- V - a - 1]J \P1P2\x + v + a][x + v - a - l][x - v - a - 1] 



At this point we can follow the Bethe ansatz procedure exactly as before. This yields the eigenvalue 
expression for the quantum transfer matrix 

A(^;) = 

C /2 A/* 



PiP2[x + v + a][x + v-a - l][x -v-a-l]J fJ-PsK " ^ + "/2 + 1] 



P3[x + v][x-v] 



PiP2Pa[x + v + Q\[x-v-a-l] 



nf P2Pi[ui - V-a/2] \y\( P^P^ ~ ^ + "^/^ + ^/2] ^ 
,=1 \P3[u^ -v + a/2 + 1]) Vp2P4 K- -v + a/2 + 1/2]) 



p4x + v][x- v] ^^^^Y\f P2K - V- q/2] \ -pt K' - + a/2 - 1/2] 

11 I U,. _ _L n,/9l i 11 



PiP2P3[x + V + a][x- V - a - 1]J [ui - v + a/2] J jj-^ \P2Pi [wj -v + a/2 + 1/2] 

. -i\-M ( [x + v - (A[x-v][x-v -"4 \^^^-prP2K-^-a/2] 



,PiP2[x + f + a] [x - + «] [x - ^ - " - 1] / 7^-7 [ui-v + a/2] 
subject to the Bethe ansatz equations 

f P3[x + Ui - a/2][x - Uj - a/2 - 1] \^^^ ^ yr PsK- - Uj - 1/2] 

V p4[x + Mi + a/2][x-Mi + a/2] / f^PAwj - Ui + 1/2] 

PaX'^ Y\ P^^^i ~ + ^/2] _ TT K' -Wk-^] 



TT /^3Fi - -r i/^J ^ _-r-r 
P3 y fj- P4 [tii -Wj-\/2] 11 [u;^- - Wfc + 1] ' 



Another expression is obtained through (|2^) . It is curious that the gauge parameters pi, p2 do not 
appear in these Bethe ansatz equations. As explained earlier, equivalent forms, which we will not make 
explicit here, can be obtained through the use of the equivalence of the expressions (|56|j58|) . 
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